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XLI 


ON THE APPLICATION OF THE METHOD OF 
QUATERNIONS TO SOME DYNAMICAL QUESTIONS 


Communicated 14 July and 21 July 1845. 


[Proc. Roy. Irish Acad. vol. m (1847), Appendix, pp. xxxvi-l.] 


The following is the substance of the communication made to the Academy by Sir William 
Hamilton, on the application of the method of Quaternions to some dynamical questions: 

The author stated that, during a visit which he had lately made to England, Sir John 
Herschel suggested to him that the internal character (if it may be so called) of the method of 
quaternions, or of vectors, as applied to algebraic geometry, that character by which it is 
independent of any foreign and arbitrary axes of coordinates, might make it useful in 
researches respecting the attractions of a system of bodies. A beginning of such a research 
had been made by Sir William Hamilton in October 1844, which went so far, but only so far, 
as the deducing of the constancy of the plane of an orbit, and the equable description of areas, 
under one common formula, namely, the following: 


dp dp 
Past to const., 


from the general expression of a central force, namely, from the equation 
dp dp 
Pa IE’ l 

which asserts merely the coaxality of the vector p and the force H z> or the existence of one 
common line along which this vector and this force are (similarly or oppositely) directed. 

Since the suggestion above acknowledged was made, Sir William Hamilton has proposed to 
himself to express by an equation, on the principles of the method of vectors, the problem of 
any number of bodies attracting according to Newton’s law: and has arrived at the formula 


da m+Am | (A) 
which may also be thus written, di?  —Aay(— Aa?) ’ 

da _ £ m' ) 

d? ~ (a-a) y- (a -a'y (B) 
and from which he has deduced anew the known laws of the centre of gravity, of areas, and of 
the vis viva, under the forms: d\2 

i (3) x.ma=0; (c) 
dt 
d da da 
Ar (D) 
d.. m/da\?, d mm’ 
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a is the vector and m the mass of one body; «’ and m’ of another; & sums for the system; t is 
the time, d the characteristic of differentiation; A (where used) is the mark of finite differencing. 

To illustrate the method of treating equations of such forms as these, let us consider briefly 
the problem of two bodies, or of one body, as it presents itself, in the method of quaternions, 
with Newton’s law of attraction, coordinates being not employed. The differential equation 
may be thus written, da M 


de a a |( — a?) Za)’ 
a being the vector of the attracted body, drawn from the attracting one; ¢ the time; d the _ 
mark of differentiation; and M the attracting mass, or the sum of the two such masses. This 
equation gives diy aie 
a —-—7a 
df df 
which expresses merely that the force is central; and gives by integration a result already 


alluded to (as independent of that function of the distance which enters into the law of 
attraction), namely, 


(1) 


= (2) 


po coc 8 
2dt di2 : 
the constant # being a new vector, perpendicular in direction to the plane of the orbit, and in 
magnitude representing the double of the areal velocity, which velocity is thus seen to be 
constant, as also is the plane. For we have at once, by (3), 


df=0; (3) 


ap+pa=0, (4) 
implying that the variable vector æ is perpendicular to the constant vector 7; and also 
[te.aa- da.a)=2f(t—t,), (5) 


if t be the value of t at the commencement of the integral. 
Make now, to distinguish between the length and direction of the vector, 


a=n, r= a) P=—1; (6) 

we shall have da=r.de+dr.t, (7) 
and because r and dr are scalar (or real) quantities, 

a.da=r%.di—r.dr, da.a=r?de.c—r.dr: (8) 

2 3 
therefore f.dt=3(a.da—da.a) = (¢.de—de.t) =r. du, (9) 
observing that the equation a 1 E E A (10) 
The fundamental equation (1) of the problem becomes, by (6) and (9), 
dda M dM 
(11) 


di dtr dt B’ 
(in which last member the order of the factors is not indifferent), and therefore gives, by inte- 
gration, since f as well as M is constant, 


ee: t.; 12 

ae = const. ; (12) 

or, as we may also write it, = =e, de=0 (13) 
$ S E ann ; 
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We have, consequently, by (6) and (4), 
da £ da fii 


as te. Wk Oo Os (14) 

and finally, by (3), ae+ex+2r=2p, (15) 
esq 

if pak, (16) 


the constant p being here not only a scalar but an essentially positive quantity, because the 
force is supposed to be attractive, or M>0, while £?<0. The equation (15) thus obtained, 
contains the law of elliptic, parabolic, or hyperbolic motion. For if we make (by way of 
comparison with known results), Me) =e, (17) 


and (a, —e)=2v, (18) 


(æ, —€) denoting here the angle between the directions of æ and e, we have (by the formula (a) of 


the abstract of last November), ae + 6a = 2er cose: (19) 
and therefore, by (15), r= eee , (20) 


which is the known equation of a conic section referred to a focus. The Greek letters, through- 
out, represent vectors: and the Italics, scalar quantities. 

Supposing that we had no previous knowledge of the properties of cosines or of conics, we 
might have proceeded thus to investigate the nature of the locus represented by the equation 
(15). This locus is a surface of revolution round the line ¢; because the differential of its equation 


being da.e+¢e.da+2dr=0 (21) 


if we cut it by a series of concentric spheres round the origin of vectors, the sections are con- 
tained in a series of planes perpendicular to e€; since 


dr=0, (22) 

which is the differential equation of the first series, gives, by (21), 
dae+eda=0, (23) 
which is the differential equation of the second series. To study more closely this surface of 
revolution (15), make a=y+a’, (24) 


y being an arbitrary constant, and a’ a variable vector; and since it must evidently give simpler 
and more symmetric results to suppose the vector y coaxial with e, than to make the contrary 
supposition, since we shall thus place the origin of the new vectors «’ upon the axis of revolution 


of the surface, let S ey—ye=0, or y=ge, (25) 


g being an arbitrary scalar, to be disposed of according to convenience. Equations (24) and (25), 
combined with (6) and (17), will give, for every point of space, 


— a= —(a—y)P =r? + ge" + glac + ea); (26) 
and therefore, for every point of the locus (15), 
— a? =r? — 2gr + ge? + 2gp. (27) 
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The second member of this last equation may be made an exact square, by assuming 


ge? + 2gp=g7, thatis, g= i P= 2a; (28) 


Pp 
l— e? 


the scalar quotient =a, or the transformation p=a(1—e?), (29) 


being thus suggested to our attention; and with this value of g we shall have, by (27), 
—a'%*=(2a—7r)?, (30) 
that is, é 2a = ,/(— a?) + ./(—a’2); (31) 
so that either the sum or the difference of the distances of any point of the locus (15) from the 
two foci of which the vectors are respectively 0 and 2ae, is equal to the constant 2a. It is not 
difficult to prove that the upper or the lower sign is to be taken, in the formula (31), according 


as e?is < or > 1. For the case e?= 1, the recent transformation fails. 
Again, to find whether the locus has a centre, we may make 


a=y'+d=9'e+, l (32) 


g' being a new disposable scalar, and 6 a new variable vector; and, after having cleared the 
equation (15) of the radical r or 4/(— «?), by writing it as follows, 


2 
a+ (p- EF) mill (33) 
2 
we get 0=('e+8)+(p+g'e—- 5") 
2 
=+ (=") +g” (de +e) +(p+g’e?)?—g’2e?; (34) 
if we make for abridgment g’ =9 —(p+g’e?). (35) 


If y’ or g'e is to be the constant vector of the centre of the locus, it is necessary that to every 
variable vector, ô, which satisfies the equation (34), should correspond another vector — ô, 
equal in length but opposite in direction, and satisfying the same equation; therefore the terms 
g" (de +e) must disappear, and we must have 


g'=0, g'=-*4=4, (36) 


the constant a being thus suggested by the search after a centre, as well as by the search after 
a second focus. Making then g’ =a in (34), we find the following equation of the surface, when 
referred to its centre, 


2 
0=8?+ (=) +ap; (37) 
in which ap =a?(1—e?) =a*(1+e?). (38) 
And because in general, for any two vectors ô, e, the following relation holds good, 
2 ie e 
("=") - (~ 3 2) + 5%, (39) 


we may write the equation (37) under the form 


0=(146) (+a?) + ( (40) 


de— e) 
3 à 
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This last equation shows that 
when de—ed=0, then 62+a?=0; (41) 


that is to say, when ô is coaxial with, or parallel to e, or, in other words, when the vector from 
the centre coincides (in either direction) with the axis of revolution of the surface, its length is 
= +a, according asais > or < 0. 

The equation (37) shows that 


when de+ed=0, then 6%+a7(1+e6?)=0; (42) 


if therefore e? be > —1, that is, if e?< 1, the length of every vector drawn from the centre 
perpendicularly to the axis of revolution will be 

\(—8%) =a(1—e2) =, (43) 
b being a new scalar quantity; but if e?> 1, e?< — 1, 1+e?<0, then we shall have, by (42), the 
absurd result of a VECTOR 6 appearing to have a POSITIVE SQUARE: whereas it is a first principle 
of the present method of calculation, that the square of every vector is to be regarded as a negative 
number: which symbolical contradiction indicates the GEOMETRICAL IMPOSSIBILITY of drawing 
from the centre to any point of the locus, a straight line which shall be perpendicular to the axis of 
revolution, in the case where e2>1. The locus has, in this case, two infinite branches enclosed 
within the two branches of the asymptotic cone which has for its equation 


2 
and nowhere penetrates within that inscribed spheric surface, which has for its equation 
d+a2=0, (45) 


though it touches this last surface at the two points where it meets the axis of revolution. On 
the other hand, when e?< 1, the locus is entirely contained within the spheric surface (45), 
touching it, however, in like manner in two points upon the axis of revolution. A finite surface 
of revolution (the ellipsoid) might thus have been discovered, of which each point has a con- 
stant swm of distances from two fixed foci; and an infinite surface (the hyperboloid), with two 
separate sheets, of which each point has a constant difference of distances from two such foci: 
and all the other properties of these two surfaces of revolution might have been found, and 
may be proved anew, by pursuing this sort of analysis. A third distinct surface of the same 
class, but infinite in one direction only (the paraboloid), might have been suggested by the 


observation that the reduction to a centre fails in the case e?=1, e?= —1. Its equation may 

be put under the form (ea" — x"e)? = 4p(ea" +a"e), (46) 
. n E > 

by making “=a -*, (47) 


so that a” is the vector from the vertex: and it lies entirely on one side of the plane which 
touches it at the vertex, namely, the plane 


ca” +a"e=0. (48) 


In general whatever e or e may be, and therefore for all the three surfaces, the length of the 
focal vector perpendicular to the axis is p; for, by (33), if we make 


cea +aec=0, (49) 
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we get a?+ p*?=0,. (50) 
Indeed (15) then gives r= p. 
Since a+r2=0, a.da+da.a+2rdr=0, (51) 


the differential equation (21) of the locus (15) may be put under the form 
(re—a)da+da(re—a) =0; (52) 


thus shewing that the vector re—« is perpendicular to the differential d of the focal vector «, 
or that it is parallel to the normal to the locus, at the extremity of that focal vector. That 
normal, therefore, intersects the axis of revolution in a point, of which the focal vector is re; 
the position of the normal is, therefore, entirely known, and every thing that depends upon it 
may be found, for the particular surfaces of revolution which have been here considered. For 
example, in the ellipsoid, the vector of the second focus, drawn from the first, has been seen 
to be 2ae; if, then, we make Qa—r=r', (53) 


so that r’ denotes the length of the second focal vector, drawn to the same point as the first 
focal vector, of which the length is r, we have — r'e for the second focal vector of the inter- 
section of the normal with the axis; the normal, therefore, cuts (internally) the interval between 
the two foci, into segments proportional to the two conterminous focal distances of the point 
upon the ellipsoid, and consequently bisects the angle between those focal distances. Again, 
if we divide the expression re — a by the scalar quantity r, and multiply the quotient by a, we 
find that e—1 and ae—a are also expressions for vectors in the normal direction; and because 
ae is the focal vector of the centre, while — az is a radius of the circumscribed sphere, opposite 
in direction to the focal vector of the point upon the ellipsoid, we see that if the focal vector of 
the extremity of this radius of the sphere be prolonged through the focus, it will cut perpendicu- 
larly the tangent plane to the ellipsoid. Again, the expression 


T=A(€+l)—a=(a—r)t+ae, (54) 


is easily seen to denote here a vector perpendicular to æ — re, and therefore to the normal, 
because aT +Ta=r(eT +7¢)=2(ap—rr’); (55) 
but 7 is also in the same plane with « and e, and therefore is a vector parallel to the tangent to 
the elliptic section of the locus made by a plane passing through the axis of revolution; az is 
therefore the central vector of a point upon this tangent, because a — ae is the central vector 
of the point of contact; and the central vector of the second focus being ae, we have ai—ae as 
an expression for the second focal vector of the same point upon the tangent; this second focal 
vector is therefore parallel to the normal, because ¿— e€ is parallel thereto, and, consequently, 
it is the perpendicular let fall from the second focus on the tangent line or plane: and the foot 
of this perpendicular is thus seen to be at the extremity of that radius of the circumscribed circle 
or sphere, which is drawn in a direction similar (and not, as lately, opposite) to the direction 
of the first focal vector of the point on the ellipse or ellipsoid. We see, at the same time, that — T 
is a symbol for the projection of the second focal vector upon the tangent line or plane; from 
which we may infer, by (55), that the product of the lengths of the two projections of the two 
focal vectors on the tangent is =7r’ —ap, and therefore that it is less than the product rr’ of 
the lengths of those two vectors by the constant quantity ap, or b?, which constant must thus 
be equal to the product of the lengths of the projections of the same two vectors on the normal, 
so that we may write the equation PP’ =ap=.?, (56) 
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if P and P’ denote the lengths of the perpendiculars let fall from the two foci on the tangent, 
while b is the axis minor of the ellipse. Analogous reasoning may be applied to the hyperbola, 
or to the surface formed by its revolution rounds its transverse axis. Most of the foregoing 
geometrical results are well known, and probably all of them are so: but it may be considered 
worth while to have briefly indicated the manner in which they reproduce themselves in these 
new processes of calculation. 
The vector drawn from the focus first considered to any arbitrary point upon the normal, 
may be represented by the expression 
v=(l—n)a+nre, (57) 


in which v is an arbitrary scalar; and if this normal intersect another normal infinitely near it, 
then we may write, as the expression of this relation, 


0=dv=(1—n)da+nedr+ (re—«a) dn: (58) 


comparing which differential equation with the forms (52) and (21) of the differential equation 
of the surface of revolution (15), we can eliminate the scalar differential dn, and deduce for n 


itself the expression da? 
S Ee: VARA (59) 
dr? + da?’ 
One way of satisfying these conditions is to suppose 
n=1, dr, v=re; (60) 


which comes to considering the intersection of the given normal with the axis, and therefore 
with the other normals from points of the same generating circle of the surface of revolution: 
and this intersection is accordingly one centre of curvature of that surface. The only other 
way of obtaining an intersection of two normals infinitely near, is to suppose, by (58), the 
element da coplanar with « and e, or to pass to consecutive normals contained in the same 
plane drawn through the axis; that is to say, the other centre of curvature of the surface is the 
centre of curvature of its meridian. The length of the element of this meridian, that is the 
length of da, is denoted by the radical ,/(—d«’*), because the differential da is a vector; and 
the length of the projection of this element on the focal vector is + dr = ./(+dr*), because dr 
is a scalar differential: therefore the length of the projection of the same element on a line 
perpendicular to the focal vector, and drawn in the plane through the axis, is denoted by this 
other radical, ,/(—da?—dr?); but the length of this last projection is evidently to the length of 
the element itself, as the length P of the perpendicular let fall from the focus on the tangent is 
‘to the length r of the focal vector of the point of contact; such, therefore, is, by (59), the ratio 
of n- to 1, if the scalar n, in the equation of the normal (57), receive the value which corresponds 
to the centre of curvature of the meridian; therefore we have 


rt ae (61) 


N denoting the length of the portion of the normal which is comprised between the meridian 
and the axis, and R denoting the length of the radius of curvature of the meridian. The projec- 
tion of this radius on the focal vector is evidently the focal half chord of curvature, of which 
half chord the length may be here denoted by C; we see then that if we again project this half 
chord on the normal, the result is the normal itself, that is the portion N, because this double 
process of projection multiplies R twice successively by n-*; and if, once more, the normal be 
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projected on the focal vector, the third projection so obtained is equal in length to the semi- 
parameter p, because, by (15) and (16), 


U(ré—a) + (re — æ) t= 2p; (62) 
hence \(RN)=C=np= — = aie í (63) 


that is, for any conic section, the geometrical mean between the radius of curvature and the 
normal is equal to the harmonic mean between the two focal distances; of which distances the 
second, namely 7’, is to be regarded as negative for the hyperbola, and infinite for the parabola, 
and the harmonic mean determined accordingly. We have also, for every conic section (if r'a~! 


be suitably interpreted), RIO İN W 

no ETEA tal 
so that the semiparameter, the normal, the focal half chord of curvature, and the radius of 
curvature, are in continued geometrical progression: and the analysis may be verified, by 
calculating directly, on the same principles, the length of the normal, as follows: 


N = J{-(re—a)*} = /fr2(e2? +1) + 2rtp—n)}=, | (P=) ; (65) 


The general relation of a conic section to a directrix is an immediate geometrical consequence 
of the equation (15), which has been here (in part) discussed, and may be regarded as its simplest 
interpretation. Some of the foregoing symbolical results respecting such a section admit of 
dynamical interpretations also; and, in particular, the expression ai — ae, which has been seen 
to represent, both in length and in direction, the perpendicular let fall from the second focus 
on the tangent, may suggest, by its composition, what is, however, a more immediate con- 
sequence of the equation (12), that in the undisturbed motion of a planet or comet about the 
sun, the whole varying tangential velocity may be decomposed into two partial velocities, of which 
both are constant in magnitude, while one of them is constant in direction also. The component 
velocity, which is constant in magnitude, but not in direction, is always in the plane of the 
orbit, and is perpendicular to the heliocentric radius vector of the body; the other component, 
which is constant in both magnitude and direction is parallel to the velocity at perihelion; and 
the magnitude of this fixed component is to the magnitude of the revolving one in the ratio of 
the eccentricity e to unity. The author supposes that this theorem respecting a decomposition 
of the velocity in an eccentric orbit is known, though he does not remember having met with 
it; but conceived that it might properly be mentioned here, as being a very easy and immediate 
consequence of the present analysis: respecting the general principles of which analysis, the 
reader is requested to consult the Abstract, already referred to, of the communication* 
of November 1844, printed at the commencement of the present volume of the Proceedings of 
the Academy. There is no difficulty in deducing, on the same principles, from formulae of the 
present paper, the known differential equation, 


Gen, 7 


which connects the radial component of velocity with the heliocentric distance and the time, 
and may be integrated by the usual processes. 


(64) 


* [See XVI.] 
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